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ABSTRACT 

Recent models of carbon ignition on accreting neutron stars predict superburst ignition depths that are an 
order of magnitude larger than observed. We explore a possible solution to this problem, that the compact stars 
in low mass X-ray binaries that have shown superbursts are in fact strange stars with a crust of normal matter. 
We calculate the properties of superbursts on strange stars, and the resulting constraints on the properties of 
strange quark matter. We show that the observed ignition conditions exclude fast neutrino emission in the 
quark core, for example by the direct Urea process, which implies that strange quark matter at stellar densities 
should be in a color superconducting state. For slow neutrino emission in the quark matter core, we find that 
reproducing superburst properties requires a definite relation between three poorly constrained properties of 
strange quark matter: its thermal conductivity (K), its slow neutrino emissivity (e„ ~ Q v x Tg) and the energy 
released by converting a nucleon into strange quark matter (2sqm)- 

Subject headings: accretion, accretion disks-X-rays:bursts-stars:neutron 



1. INTRODUCTION 

Strange stars are compact stars that are almost entirely 
made of deconfined quark matter, long ago suggested as an 
alternative to neutron stars following the proposal that strange 
quark matter (SQM) is the groun d state of hadronic matter 
( Witter] fl984t see also Iltohl fl970) (For a review see Weber 
2005). Strange quark matter with a density of ~ 5 x 10 14 
g cm -3 could exist up to the surface of a strange star, in 
which case the star wo uld look, and evolv e, very differ- 
ent from a neutron star (Page&Usov 2002). However, a 
stra nge star can also have a crus t of normal baryonic mat- 
ter ( Alcock. Farhi & Olinto 1986), in which case its surface 
emission properties would be the same as those of a neu- 
tron star. The bulk properties of both kinds of stars, such 
as radius and moment of inertia, are also similar in the ob- 
served mass range 1 <M/Mq < 2, and so in general it is dif- 
ficult to distinguish between s trange stars and neutron stars 
jHaensel. Zdunik & Schaefferll986l) . To date, there is no real 
evidence against or in favor of the strange star hypothesis 
dPrakash et alJl2003l) . mostly due to the fact that many prop- 
erties of quark matter are poorly understood, allowing great 
freedom to vary parameters of the models to reproduce ob- 
served properties. An observable, and observed, characteristic 
of compact stars that would allow us to distinguish between 
these two very different kinds of objects has not yet been iden- 
tified. 

Superbursts are rare, extremely energetic, and long duration 
Type I X- ray bursts o bserved from low mass X-ray binaries 
(LMXBs) ( Kuulkers 2004). They have been modeled as being 
due to unstable thermonuclear b urning of carbon on the sur- 
face of an accreting neutron star (gumming & Bildste rl2001t 
IStrohmaver & Brownll20 02 ) . The accreted hydrogen and he- 
lium burns shortly after its arrival on the stellar surface, pro- 
ducing a mixture of heavy elements and carbon (X c ~ 10% 
carbon by mass) that accumulates an d later ignites. Brown 
(2004) and Cooper &~Naravanl ( 120051) showed that the igni- 
tion conditions for superbursts are sensitive to the thermal 
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properties of the stellar interior, in particular the composition 
of the crust, and the neutrino emission from the core. This 
opens up a new way to probe the interiors of the compact stars 
accreting in LMXBs. 

However, calculations so far suggest that carbon ig- 
nition models, on a neutron star, have difficulty in 
achieving the low ignition column depth s inferre d from 
observations {C umming & Bildsten 2001: iBrownl 2004; 
Cooper & Naravan 2005; Cumming et al. 2005). In particu- 
lar, gumming et al. (2005) carefully constrained the ignition 
depth by fitting models of the cooling of the burnin g layer 
following a superburst ( Cumming & Macbeth 2004) to ob- 
served superburst lightcurves. They then showed that the in- 
ferred ignition column depths of y = (0.5-3) x 10 12 g cm" 2 
were an order of magnitude smaller than carbon ignition 
models predict. In particular, they showed that neutrino 
cooling i n the crust from Cooper pairing of the superfluid 
neutrons (Flowers. Ruderman. & Sutherland 1976) limits the 
crust temperature to < 5 x 10 8 K, making it impossible to 
achieve ignition at the inferred depths. 

In this paper, we investigate a possible solution to this puz- 
zle, which is that the compact stars in LMXBs that have 
shown superbursts are in fact strange stars with a crust of nor- 
mal matter. These stars do not have an inner crust, and so 
neutron Cooper pair neutrino emission is absent; in addition, 
the energy release at the base of the crust from conversion 
of accreted nucleons into strange quark matter is potentially 
much greater than the energy release in neutron star crusts 
from pycnonuclear reactions and electron captures. We there- 
fore expect the accreted crust of a strange star to be hotter than 
for a neutron star, making carbon ignition at column depths 
< 10 12 g cm" 2 possible. We first describe our strange star 
models in § [2] and then present calculations of carbon igni- 
tion in §13 We conclude in §0 

2. STRANGE STAR MODELS 

We first discuss our models of strange stars with a thin 
baryonic crust, and the physical parameters relevant for car- 
bon ignition. The accumulating carbon layer is heated by the 
energy released at the base of the crust from conversion of 
nucleons into strange quark matter. Most of this energy is 
conducted into the core and escapes as neutrinos, but a small 
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FIG. 1 . — Masses of strange stars vs central density (left) and radius (right) 
for four sets of values of the coupling constant etc, bag constant B 1 / 4 (in 
MeV), and strange quark mass m$ (in MeV/c 2 ) covering a broad range of 
possible values. 
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FIG. 2. — Density profiles of the outer layers of four strange stars with 
crusts. The models have different values of m s but the same values of a c — 0.3 
and B 1 / 4 = 140 MeV. The electron chemical potential at the surface of the 
quark matter is indicated, as well as the infe rred range of carbon igni tion 
density for superbursts. The crust EOS is from Haensel & Zdunik 1 1990). 



amount flows outwards through the crust and sets the igni- 
tion conditions. The important ingredients are (1) the crust 
thickness, (2) the energy released per nucleon by conversion 
into strange quark matter, Qsqm, (3) the neutrino emissivity 
e v and (4) thermal conductivity K of quark matter. Given our 
poor understanding of quark matter these four ingredients are 
essentially free parameters and we will treat them as such, 
calculating the values required to match observed superburst 
properties. 

To obtain estimates of the possible range of values for the 
crust thickness and Qsom we consi der the strange quark mat- 
ter models of Farhi & Jaffe ( 1984). These MIT bag inspired 
models use three parameters, the bag constant B which is 
the energy density cost for having deconfined quark matter, 
the strange quark mass m s , and the color interaction coupling 
constant a c . For a broad range of values of these parameters 
strange quark matter is self bound. Figuref^shows examples 
of mass-radius and mass-central density relationships for a 
broad range of parameters. For definiteness, in our numerical 
models we choose a c = 0.3, fi=(140 MeV) 4 , and M = 1 .4 M Q . 

The presence of a crust of baryonic matter is possible be- 
cause of the positive charge of strange quark matter which 
hence must contain electrons 3 . The electrons are bound to 
the quark matter by the long range Colomb force, and leak 
out of the quark matter creating a dipolar layer of thickness a 
few thousand fms. The outward pointing electric field repels 
ions, el ectro statically suspendi ng the crust above the quark 
matter ( Alcock. Farhi & Olinto 1986). The existence of the 
supporting electric field requires that the electron chemical 
potential /i e in the crust be lower than in the quark matter. 
This sets the thickness of the crust. Accretion drives contin- 
uous energy release, because as the density at the base of the 
crust increases, the value of n e there grows and approaches its 
value at the surface of the quark matter resulting in a reduc- 
tion, or even inversion, of the electric field: ions plunge into 
the quark matter, releasing an energy Qsqm per baryon. The 
thickest possible crust, which occurs when fi e in quark mat- 
ter is larger than ~ 30 MeV, reaches neutron drip at its base: 
the dripped neutrons fall directly into t he quark matter, also 
producing an energy Gsqm- Within the lFarhi & .Taffd if 1984) 
model, for fixed B and a c , the value of [i, e at the surface of the 
quark matter core is determined by the values of m s , a larger 

3 We do not consider models having negatively charged quark matter 
which would contain positrons instead of electrons. 



m s resulting in a larger fi e . Figure |2] gives the density profile 
of the outer layers of the four strange star models we will use. 
The thickness of the nuclear crust ranges from rj 100 m to a 
maximum extent, at neutron drip, of w 400 m. 

The energy released per nucleon by converting it into 
strange qu ark matter, (?som, i s our second free parameter. 
Within the Farhi & Jaffe 1 1984) model it can take values from 
a few MeV up to almost 100 MeV, depending on the precise 
values of the model parameters (B, a c , m s ). For our numeri- 
cal simulations we vary Qsqm over a broad range 4 , from < 1 
MeV to > 100 MeV. 

The neutrino emissivity of quark matter can span the broad- 
est possible range. If quarks form a normal Fermi liquid, the 
direct Urea process, u + e~ — > d + v e with d — > u + e~ + V e (and 
the same process with an s replacing the d), has an enormous 
emissivity of the order of e^ Urca ~ 10 25 T 9 6 erg cm" 3 s" 1 , where 
Tg is the temperature T in units of 10 9 K. However, in the very 
likely situation t hat some flavo rs/colors of quarks form a color 
superconductor I Alford 2001), the neutrino emission may be 
determined by b remsstrahlung from unp aired quarks, or from 
electrons (e.g. JPage. Geppert & Webea2 005). We write these 
slow neutrino emission processes as 

e v ~ Q v x T 9 & erg cm" 3 s" 1 (1) 

and consider Q v as our third free parameter with values from 
about 10 18 up to an optimistic maximum of ~ 10 22 . 

Finally, the thermal conductivity K of quark matter is our 
last free parameter. We will consider values from K = 10 19 to 
10 22 erg cm " 1 K" 1 , which cove r the whole range of possible 
values (e.g. JPage & Uso\l2002l) . 

3. CARBON IGNITION IN STRANGE STARS 

We now present our numerical calculations of the thermal 
state of an accreting strange star, and investigate the range of 
values of parameters that result in successful superbursts, i.e., 

4 Strictly speaking, varying Qsqm means varying the EOS while we use 
only four cases of EOS and vary Qsqm for each. However our results are only 
weakly dependent on the mass and size of the star and adjusting the EOS for 
each v alue of Qsqm would not affect them. However, within the Far hi & Jaffd 
i 1984) models large values of m s result in a reduction of the maximum possi- 
ble value of Qsqm ar, d a l so an increase in fi e , i.e., thicker crusts: cases which 
allow ii e > 30 MeV, and hence a maximally thick crust, limit Qsqm to De less 
than 40 MeV. Considering that this is only an approximate model we choose 
to explore the broadest possible range of values of Qsqm but one must keep 
in mind this possible restriction of Qsqm m analyzing our final results. 
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FIG. 3. — Temperature profiles for strange stars with a thick crust reaching 
the neutron drip density, accreting at m = 0.3mEdd- We take a slow neutrino 
emissivity (eq.0 with Q v = 10 21 , but different values of K and 2sqm as 
indicated. The carbon ignition curve (with Xq = 0.2) is shown as a dashed 
line. The two vertical lines indicate the observed range of ignition depths. 
The dotted curve shows the melting line. Note that most of the crust is liquid 
(in fact, we find that thinner crusts can be completely liquid). 



carbon ignition at column depths of y w 10 12 g cm" 2 . We as- 
sume a constant accretion rate m (in nucleons cm" 2 s" 1 ) and 
follow the time evolution of the star until it has reached a sta- 
tionary state. Typical values of m are about 0.1-0.3 ;«Edd (in 
terms of the Eddington rate \-KR 2 m u m^i — 10 18 g s" 1 , R being 
the star's radius and m u the atomic mass unit). Heat is injected 
at the top of the strange quark matter core at a rate given by 
H = <2sqm x 4-TrR 2 m. Our code solves the heat transport and 
energy balance equations in full general relativistic form with 
spherical symmetry. 

A first and firm result of our simulations is that it is impos- 
sible for the crust to reach temperatures close to the required 
w 6 x 10 8 K if a direct Urea process operates in the core, 
even if suppressed down to e u ~ 10 22 r 9 6 erg cm" 3 s" 1 . There- 
fore, the neutrino emissivity of strange quark matter must be 
restricted to the slow processes of equation Q. Illustrative 
examples of temperature profiles are shown in Figure [3] as- 
suming Q v = 10 21 erg cm" 3 s" 1 , but different values of K and 
Qsqm- The model with K = 10 20 cgs and Q SQM = 100 MeV 
is too hot, and ignites at a column depth smaller than ob- 
served, whereas the model with K = 10 22 cgs and Qsqm = 20 
MeV is too cold, and ignites at a larger column than observed. 
However, we obtain ignition at the appropriate depth for ei- 
ther K = 10 22 cgs and g SQM = 100 MeV, or K = 10 20 cgs and 
£soM = 20MeV 

We can understand the temperature profiles in Figure [3] 
with a simple analytic model. Because the crust is ap- 
proximately isothermal below the ignition depth, the tem- 
perature at the crust-core interface is a reasonable esti- 
mate of the ignition temperature. In addition, most of 
the heat released at this interface flows into the core 
JCumming et al.l 120051 showed that an outwards flux of 
only Qb ~ 0.25 MeV per nucleon (m/0.3?T2Edd) _I -C Qsqm is 
needed to achieve ignition at y w 10 12 g cm" 2 ). When the 
thermal conductivity is large, this heat is rapidly distributed 
over the core, which is then isothermal at a temperature given 
by balancing the heating with the energy loss due to neutri- 
nos: 4irR 2 mQsQM = (4vt/? 3 /2>)Q v T^. We therefore find that 
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FIG. 4. — The energy release Qsqm required for carbon ignition at 
10 12 g cm -2 as a function of Q v and K. The thickness of each shaded stripe 
corresponds to the range of crust thickness in our models (Fig. [2] a thicker 
crust requires a larger Qsqm)- The criterion for ignition is T = 6 X 10 s K at 
y = 10 12 g cirr 2 . To account for the allowed range of y shown by the vertical 
solid lines in Figure l3l the upper limit of each strip can be increased by a 
factor two (corresponding to y ~ 0.5 X 10 12 g cm -2 ), or the lower limit de- 
creased by a factor two (y ~ 3 X 10 12 g cm -2 ). The dashed lines show the 
analytic estimates for the same values of K. 

T 9 = OmQsQu/RQvy^, or 

T 9 = 0.54 Q-Jif ( ( -*—) 'Vv* (2 ) 
» WJ "- 21 v 100MeVy \03thim) 6 

where Q v<2 i = G^/10 21 erg cm" 3 s" 1 and R 6 = R/10 km. This 
estimate agrees well with Figurel^for models with an isother- 
mal core (K = 10 22 cgs). 

However, when the thermal conductivity is low, the core is 
no longer isothermal, and we must estimate the temperature 
profile in the core. To do this, we use an approximate plane 
parallel model of the core. The temperature profile is given 
by pdF /dy = —e u , where F = -pKdT jdy is the inwards di- 
rected heat flux, p being the density. Assuming p, Q v and K 
are constant in the core gives d 2 T /dy 2 = Q V T^ / p 2 K. The so- 
lution is T oc (y + b)~ 2 l' ! ', where b is a constant determined by 
the requirement that all of the heat flux entering the core is 
emitted as neutrinos: wiQsoM = J £udy. Assuming the surface 
contribution dominates the integral, we find the interesting re- 
sult (mgsoM) 2 = (2/9)(RQ„T & )(KT/R); the heat flux entering 
the core is the geometric mean of the neutrino emission at the 
core surface and the heat flux transported deep into the core. 
The temperature at the core surface is 

^^'2i«"(^) 2 "(^) 2 ",(3) 

where K20 = K/10 20 cgs, which agrees well with the non- 
isothermal core profiles in Figure [3] Comparing equations 
(0 and gives the thermal conductivity of the core at 
which the temperature profile becomes isothermal as K crlt = 

4.3 x 10 21 cgs e^VeWlOOMeV^V/OJ m Edd ) 7 / 8 . 

Figure |4] shows the combinations of 2sqm, Q u , and K 
that are needed for ignition at y = 10 12 g cm" 2 . We find 
the surprising result that even values of (2sqm as large as 
100 MeV (compared to the « 1.5 MeV released by pyc- 
nonuclear reactions in the neutron star case) can give the 
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required crust temperature when K and Q v are sufficiently 
large. The shape of the curves is determined by the fact 
that, at a fixed value of Q v , as the conductivity K is in- 
creased eventually the core becomes isothermal, and the solu- 
tion then becomes independent of K. For an isothermal core, 
fixing the temperature at the boundary to be Tg = 0.7 (which 
gives ignition at the right depth; Fig- 0» in equation (|2ji gives 
GsQM.iso = 80 MeV g„,2ofl6('«/0.3 (n Edd ) _1 . Similarly, set- 
ting Tg = 0.7 in equation for the non-isothermal case gives 

GsQM.noniso = 38 MeV Q^K^ 2 (m / 3 WEddT 1 . These an- 
alytic estimates are plotted as dotted lines in Figure |4] and 
nicely reproduce the numerical results. 

4. DISCUSSION AND CONCLUSIONS 

In this paper, we have explored the scenario that superbursts 
originate on strange stars with a thin baryonic crust. Given our 
poor understanding of quark matter, it is perhaps not surpris- 
ing that we have been able to find many choices of parameters 
for accreting strange stars that reproduce the observed proper- 
ties of superbursts. Nevertheless, the models provide signifi- 
cant constraints on the properties of strange quark matter. In 
particular, we cannot reproduce the observed ignition condi- 
tions if the neutrino emission in the quark matter core is by a 
fast process such as direct Urea. The fact that direct processes 
must be suppressed implies that strange quark m atter at stel- 
lar densities should be a color superconductor ( Alfordl l2b01t 
iPa ge. Geppert & Weberl 12005). For slow neutrino emission 
in the quark matter core, we find that reproducing superburst 
properties requires a definite relation between three poorly 
constrained properties of strange quark matter: its thermal 
conductivity (K), its slow neutrino emissivity (e„ ~ Q u x Tg) 
and the energy released by converting a nucleon into quark 
matter (Qsqm)- The phenomenlogically required correlations 
are shown in Figure|4] 

If strange stars exist, it is possible that either all com- 
pact stars are strange stars, or only the most massive neutron 
stars are strange stars, depending on whether the unknown 
critical density at which the transition to strange matter oc- 
curs is reached in core-collapse supernovae or whether ac- 
cretion of matter is required to increase the neutron star cen- 
tral density to the transition point. Nevertheless, in both of 
these extreme cases, and assuming strange stars exist, the 
most massive compact stars must be strange stars. Recip- 
rocally, if the most massive compact stars can be proved 
not to be strange stars we would be able to conclude that 
strange stars simply do not exist at all. The LMXBs are 
good candidates for containing very massive compact stars, 
and hence possibly strange stars, because they are long-lived 
systems which undergo a significan t amount of mass transfer 
(as much as 0.2-0.7 M^: ICook Sh apiro & TeukolskyjQUl 
Ivan den Heuvel & Bitzarakill 19951) . although how much mat- 
ter accretes onto the compact object is an open question. Two 



recent mass measurements of millisecond radio pulsars with 
white dwarf companions (likely to be descendents of LMXBs) 
give very diff erent values of 2. 1 ± 0.2 M Q for the 3.4 ms PSR 
J0751+1807 jNice et al.|l2005l) . confirming the formation of 
very massive compact stars in LMXBs, and 1.438 ±0.024 
M Q dJacobi et alJl2005l) . for the 2.95 ms period PSR J1909- 
3744, showing that not all LMXBs will result in very massive 
stars. The high mass of PSR J0751+1807 has the very impor- 
tant consequence that if strange stars exist the EOS must be 
very stiff and, as shown in Figure^ this implies that the cen- 
tral density of strange stars does not change by much more 
than a factor two when the mass increases from 1.1 to 2.1 
Mq. In other words, all strange stars must be very similar to 
each other, and behave similarly independently of their mass, 
in contradistinction to neutron stars where various phases of 
exotic matter may appear at different densities, and hence dif- 
ferent masses, and result in very distinct behaviors. So far, 
we cannot claim that superbursts provide concrete evidence 
in favor of the existence of strange stars, although superbursts 
from strange stars do offer a "natural" explanation for the ob- 
served low ignition column depths for carbon. More theo- 
retical work on neutron star models for superbursts is needed 
before we can say definitively that a strange star solution is 
required. 

Several interesting questions remain to be addressed in fu- 
ture work on accreting strange stars. Our models assume 
steady-state accretion, but in fact several objects are under- 
going time-dependent accretion. Modeling these transient 
systems will give us extra, and very tight, constraints, par- 
ticularly on the star's specific heat. An important example 
is KS 173 1-260 which was observed to accrete for 13 years 
before transitioning into quiescence in early 2001. Its qui- 
escent flux has been measured and is too low to explain the 
super burst that was observed from this source in early 2001 
jRutledee et alJl2002t |Browr3[2004). In cases for which the 
nuclear crust is thin, 100-200 meters, we find that the whole 
crust is liquid: a nuclear ocean, electrostatically suspended 
over quark matter whose density is at least four orders of mag- 
nitude higher, must have a very distinct hydrodynamic behav- 
ior which should be investigated. Finally, finely constraining 
the properties of strange quark matter through modeling of 
these phenomena, we will be able to study the cooling of iso- 
lated strange stars. This will tell us if average compact stars, 
with masses between 1. 1-1.4 Mq, are or are not strange stars. 
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